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Abstract Efforts are made to understand the timing behaviors (e.g., the jumps 
of the projected pulsar semimajor axis in the periastron passages) observed via 
13-year timing for PSR B1259-63. In the first, planet-like objects are suggested 
to orbit around the Be star, which may gravitationally perturb the (probably low 
mass) pulsar when it passes through periastron. Second, an accretion disk should 
exist outside the pulsar's light cylinder, which creates a spindown torque on the 
pulsar due to the propeller effect. The observed negative braking index and the 
discrepant timing residuals close to periastron could be relevant to the existence 
of the disk with varying accretion rates. Third, a speculation is presented that the 
accretion rate may increase on a long timescale in order to explain the negative 
braking index. 

Key words: pulsars: individual: PSR B1259— 63 — binaries: general — stars: 
neutron — stars: early-type 



1 INTRODUCTION 

Pulsars are clocks with great precision, with which one can probe astrophysical processes. The 
pulsar clock technique has been very successful in studies of dynamical systems (both general 
relativistic and Newtonian) and pulsar interiors (glitches) . Thirteen years of timing data have 
been obtained for the unique pulsar PSR B1259-63, which has a massive companion with a 
circumsteilar disk (Wang et al. 2004). Unfortunately, as concluded by Wang et al. (2004), a 
plausible mechanism to account for the observed timing behavior still has not been offered. 

PSR B1259-63 was discovered in a survey of the Galactic plane by Johnston et al. (1992), but 
determining its timing properties has not been easy. Manchester et al. (1995) fitted the timing 
data by introducing AP/P ~ 10~ 9 at each periastron, which they attributed to propeller torque 
spindown. Wex et al. (1998) presented a timing model with the inclusion of changing to and 
x = d/di[a p sini] terms, which still failed to accurately describe the 13-year dataset (Wang et 
al 2004). Jumps in projected semimajor axis Aa; = A(a p sini) at periastron were introduced 
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by Wang et al. (2004). The rms residual in the orbit-jump model is only about half that in the 
i/-i>-jump model, even though the free parameters in the later model are four more than in the 
former one. 

PSR B1259-63 may accrete from the dense circumstellar environment when it passes though 
periastron. Observations of the pulsar enable one to probe the circumstellar medium of an 
accreting pulsar. Recent studies show that propeller disks around such compact stars may play 
an important role for various astrophysical objects, including anomalous X-ray pulsars and soft 
7-ray repeaters (AXP/SGRs, e.g., Chatterjee et al. 2000) and compact center objects (CCOs, 
Xu et al. 2003). Recently, evidence for a fossil disk around an anomalous X-ray pulsars, 4U 
0142+61, has been proposed (Wang et al. 2006) through detecting emission at two mid- infrared 
bands (4.5/zm and 8.0/xm) 1 . PSR B1259-63, as an ideal "laboratory", can be studied in great 
detail, and is the only known case of a pulsar's disk with material replenishment (from the 
wind-disc of the Be star). 



2 THE I-JUMP SCENARIOS PROPOSED 

The inclination of the system to our line of sight angle i = sin - {[/(M p )(M p + M c ) 2 ] 1 / 3 /M c } is 
~ 35° (i = 32° ~ 36° for M p /M Q = ~ 1.4) for the companion Be star with mass M c = 1OM 
and the mass function f(M p ) = 1.53M Q . The total orbital angular momentum L for a Keplerian 
orbit reads 

l ~ m p mJ2£=$p = mJ Gil - e2)M ^ ~- 

P C V Mp+M c P V S1BI (1) 

- 3 x 10 53 MpiM c Y 2 g cm 2 /s, 

for the PSR B1259— 63 system when M p <C M c (i ~ 35°), where a p is the semimajor axis and 
the projected pulsar semimajor axis x p = a p smi, M p — M M p i, M c — 1QMqM cW ("p" and 
"c" denote the pulsar and the company Be star, respectively). 

The spin angular momentum S c of the Be star can be approximated as 

S c ~ 0AM c R 2 c n c ~ 0.4VG/M c 3/2 i?c /2 
~ 8 x 10 52 /o.5M c 3 / 2 i4 /2 g C m 2 /s, 



where the Be star is rotating at / times of the centrifugal breakup frequency flk ~ w 'GM C / 'i? 3 , 
/ = O.5/0.5, the Be star radius R c — §RqR C 6- It is evident that the magnitude of L and that of 
S c are comparable. It would not be difficult to change L (while keeping |L| almost the same) 
via occasional spin-orbit coupling (i.e., transferring orbital to spin angular momentum, or vice 
versa) . 

As noted by Wang et al. (2004), it is not likely that the x-jumps are changes in semimajor 
axis a p . Rather, assuming a constant a p , then a Ax = A(a p sini) requires an inclination angle 
change per periastron passage 

Ai = (Ax p /x p ) tanz - 3° x 10~ 4 Aa;pio, (3) 

where Ax p = 10ms c Ax pW , with c the speed of light. If |L| is the same, the orbital angular 
momentum change AL, which is nearly perpendicular to L, is 

AL ~ 5.4 x 10" 6 Aa;pioL - 10 48 Ax w M pl M^ g cm 2 /s, (4) 

which shows AL <C L. 



1 However, no emission is detected at and 70/im, with limits of 0.05 mjy and of 1.5 mjy, 

respectively (Bryden et al. 2006). These two observations are consistent with Xu (2005), where at least 
part of the millisecond pulsars are suggested to be of supernova origin via accretion-induced collapses. 
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No Ai can occur if the mass distribution is symmetric with respect to the orbital plane. 
Note that Ai should have only one sign as long as the asymmetric pattern of mass-distribution 
is almost the same (e.g., the case in which an asymmetric pattern originates only from an angle 
between the spin axis of Be star and the orbit normal). One must therefore find mechanisms 
that can cause an asymmetric, time-variable mass distribution. 

Planets around the Be star? Neither observational nor theoretical conclusions of massive 
star formation and evolution are certain. Nevertheless, a general view for B-spectral class stars 
could be as follows: After gravitational collapse of a molecular cloud, a massive star may form 
inside the cloud, with a dense, dusty circumstellar disk and strong bipolar outflow. The star 
could then become a so-called young Herbig Be star if it is still embedded in a nebula (Fuente 
ct al. 2003). It would then evolve to be a normal Be star (like SS 2883, the companion to PSR 
B1259-63, without a nebula). The outflow is weak (or there is a tenuous polar wind), and its 
disk becomes less and less massive. Finally, the star becomes a B star when the disk disappears 
(like the companion of PSR J0045-7319). 

One may speculate that planet-like objects (planets, proto-planets, or planetesimals) could 
form in the dusty disks of a Be star after the first supernova of binary, with highly eccentric 
orbits. The number of such objects and thus the total mass decrease with time due to their 
being captured, disintegrated, or evaporated by their hot Be stars (or via the interaction with 
the disk- wind). Evidence for planetesimal (with asteroidal size) infall onto a very young (~ 10 
years) Herbig Be star LkH Q 234 has recently been reported (Chakraborty, Ge & Mahadevan 
2004). More than 100 extra-solar planets, with masses ~ (1 — 10)Mj up it e r 5 have been discovered 
by radial- velocity surveys (Tremaine & Zakamska 2004). Alternatively, planets formed before 
the first supernova could probably be residual if the supernova occurs far way its company at 
a very highly eccentric orbits. 

There might be many planet-like objects (planets, proto-planets, planetesimals, or even 
brown dwarfs) orbit around the Be star; but the exact number could hardly obtained by ob- 
servation or calculation. PSR B1259-63 would be perturbed by the gravity of one or more such 
planet-like bodies near its periastron if they orbit the Be star with their semimajor axis being 
much smaller than the pulsar's, a p . However, if the pulsar is though to have a conventional, 
normal mass ~ IAMq, the perturbation should be negligible. Nonetheless, pulsars could be 
quark stars having masses much smaller than ~ Mq (Xu 2005). Hence the perturbation might 
be strong enough for the observed x-jumps if PSR B1259-63 is a low-mass quark star. The 
orbital angular momentum of a planet around the Be star is i p i an ct ~ 2 x 10 47 MjO^qq g cm 2 /s, 
where the planet mass is Mj x Mj up itor and its semimajor axis is aioo x lOOi?0. From Eq.Q, 
one sees AL < L p \ anct if Mj -C M p < Q.2Mja\^Ax^Q Mq, and a significant x-jump of a low- 
mass quark star is possible through the gravitational perturbation of planet-like objects. An 
encounter of a planet and the pulsar per periastron passage would produce (1) a large exchange 
of angular momentum between the orbits (thus an x-jump) and (2) a large change in the orbital 
elements of the planet (thus a higher eccentricity of the planet). The migration inwards due to 
the interaction of the planets with the disk would increase the possibility of the encounters and 
may decrease the eccentricity, but the pulsar might in turn "eject" planets after encounters. 

It may be hard to believe that planets would exist around a Be star since it is conventionally 
thought that planets form by accretion over ~ 10 7 years. Nevertheless, planet formation could 
be sped up due to fragmentation of an unstable disk (Boss 2003). Another possibility that can 
not be ruled out, is that the planets could be strange planets (i.e., strange quark matter with 
planet masses) which were born during the first supernova explosion. This kind of planet can 
not evaporate, and thus would exist for a long time. 

Asymmetric mass distribution of the Be star? Besides the suggestion above, two more sce- 
narios are speculated for x-jumps. The periastron separation between the pulsar and the Be 
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star is r p ~ 24i? c ^> R c . If the Be star has stellar oscillations leading to multipoles of mass 
distribution, then the torque acting on the pulsar may result in an £-jump per periastron pas- 
sage. As no torque is related to the monopole term, the quadrupole one would be important for 
cc-jumps. It worth noting that the oscillations could be enhanced if they are excited resonantly 
(e.g., when the ratio of the orbital period to an oscillation period is an integer) by the orbit- 
ing pulsar (Witte & Savonije 1999). Chaotic orbital dynamics near the periastron could also 
be possible in this case (Mardling 1995), which would result in a random transfer of angular 
momenta between comparable values of L and S c (En. l|ll2|l ). Additional study of the structure 
and evolution of massive stars, especially their oscillation behaviors, might further elucidate 
the periastron puzzle in this scenario. 

Asymmetric mass distribution of the disk of SS 2883? The existence of the Be star disk is 
an obvious asymmetry with respect to the orbital plane. But this can not cause sign-changed 
x-jumps if the disk is almost homogeneous. However, if a density-wave pattern can be excited 
(e.g., by the pulsar or planets) in the disk, like the case of Saturn's ring structure affected by 
the satellite Mimas, then angular momentum may transfer between the orbiting pulsar and the 
disk in the periastron passages, which might result in £-jumps. But it seems difficult to limit 
the duration of this effect (if it exists) only to the period of the pulsar passing near the Be 
star's disk. 

3 PROPELLERS IN A RADIO PULSAR PHASE? 

One strange feature of this system is that the braking index is about —37 after 13 years of timing 
of PSR B1259-63, while the indices are between 1.4 and 2.9 (Ss 3) for 6 pulsars, with great 
certainty (Xu & Qiao 2001; Livingstone et al. 2006). It is true that the data- fitting is improved 
significantly when ir-jumps at periastron passages and a pulsar glitch near MJD 50691 are 
included, whereas a few residual points close to the periastron are still discrepant systematically. 
This could be due to errors in the DM correction (Wang et al. 2004). An alternative explanation 
presented below for the unusual braking index as well as the timing residuals is that the pulsar 
may capture SS 2883's circumstellar matter to form a pulsar disk which results in a torque 
being antiparallel to the pulsar's spin (the viscosity timescale of accretion disk could be order 
of years). 

By observing the un-pulscd radio emission at four different bands, Johnston et al. (2005) 
conclude that the pulsar passes through the dense circumstellar disk of the Be star just before 
and just after periastron. It iseems intuitive that accretion would occur when the pulsar passes 
through the disk. Accretion of in-falling matter with negligible angular momentum was consid- 
ered by Kochanek (1993) and Manchester et al. (1995), in which the pulsar loses its angular 
momentum in order for the in-falling matter to obtain enough angular momentum to escape. 
Besides spinning down the pulsar, the no-angular-momentum accretion may also affect pulsar 
orbit parameters through frictional drag (Wex et al. 1998). However, if the material captured 
by the pulsar carries significant angular momentum, disk accretion around the pulsar occurs. 
To make this more precise, one should calculate the circularization radius to see if the accreted 
matter has sufficient angular momentum for disk formation. 

Accretion disk formation around PSR B1259-63. The cylindrical radius, r acc , within which 
a pulsar with mass M p can gravitationally capture material with a relative velocity i> re i can be 
estimated as (e.g., Lipunov 1992) 

r acc ~ GMp/vli ~ 1.3 x 10 12 M pl «-f 100 cm, (5) 

where Urelioo = w roi x 100 km/s. Note that the radius of PSR B1259-63's light cylinder is 
only ric = cPj (27r) = 2.3 x 10 s cm. If the accreted material is from the Be star's wind with 
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velocity v w ~ 500 km/s, the circularization radius (Frank et al. 2002) in this case is only 
R^i TC — G 3 M 3 Sl 2 rb /u 8 ~ 200M 3 X cm, which is much smaller than the pulsar's radius, where 
^orb = 27r/P or b. Therefore a disk around the pulsar is not likely to form by accretion from the 
Be star's wind 2 . 

However, if the gas in the Be star's disk moves in a Keplerian orbit and the pulsar has an 
orbital velocity near periastron (r p ~ 24i? c ~ 10 13 cm) of v pcT i a — y/GM c (2/r p — sin35°/a;p) ~ 
157 km/s, the accretion radius r acc is different for material inside and outside periastron. The 
Keplerian velocity periastron is Vk{r p ) = y/GM c /r p ~ 115 km/s, but the sonic speed of the line- 
emission disk with temperature ~ 10 4 K is only v s ~ 10 km/s <C Vk- We will then neglect the 
pressure of disk material in the following consideration. Assuming that the pulsar orbits in the 
same direction as the circumstellar medium around the Be star, for M p = lJVf©, the Keplerian 
velocity at the accretion boundary outside or inside the periastron is Vk = y GM C / (r p ± r acc ). 
Using Eq.© and w ro i — Vk — v pcr i a , one can find r acc ~ 3.8 x 10 12 cm and ~ 6.7 x 10 12 cm for 
outside and inside the periastron, respectively. The specific angular momentum 3 is then order of 
I ~ r acc v rc i, and the circularization radius is then R^ ilc — l 2 /(GM p ) ~ 10 11 cm for r acc ~ 10 12 
cm, being much larger than r m (Fig.l, r\ c ~ 10 8 cm). One can still have R^ irc ~ 10 U ^ r ic 
even for a small pulsar mass of M p = 10 _3 Mq (r acc ~ 7.7 x 10 9 cm and ~ 7.2 x 10 12 cm for 
outside and inside the periastron, respectively, in this case). Therefore an accretion disk can 
form around the pulsar if the circumstellar disk of the Be star is Keplerian, even without the 
inclusion of magnetospheric interaction 4 , since the circularization radii are usually much larger 
than the magnetospheric ones. 

Propeller torque and the pulsar's spindown. For simplicity, we assume that the angular 
momentum of the pulsar's disk is parallel to the spin axis of PSR B1259-63. If the circumstellar 
matter around the Be star is a hydrogen plasma, the mass density would be (Wex et al. 1998) 
p(r) ~ 7.5 x 10~ 12 (r/i? c )~ 4 - 2 g/cm 3 , where the Be star's radius R c ~ 6i£ . When the pulsar 
passes through the dense disk near periastron, the accretion rate could be 

Mo ~ vrl cc v, cl p{2AR c ) ~ 5 x 10 15 M 2 X g/s, (6) 

where we choose v rc i ~ y (vk — w pC ria) 2 + ^ 2 ~ 50 km/s for the following estimates. The pulsar 
could then accrete AM ~ 9 x 10 21 Mp X g for ~ 20 days near each periastron passage. An effective 
accretion rate averaged over the orbital period would then be M ~ (20 days/P or b)M . 
The magnetic momentum can be estimated from f2 = 2ir/P and Q, 

/i ~ {^f^pRD^ „ 3 x 10 29 M p 5 / 6 G cm 3 , (7) 

for a pulsar with radius R p and mass M p (averaged density p = 3M p /(47ri? 3 ) ~4x 10 14 g/cm 3 
for quark stars with mass <^ Mq). The magnetospheric radius would then be 

r m ~ (2GM p )- 1 / 7 M 4/7 M- 2 / 7 ~ 1.2 x 10 13 M^ /3 M" 2 / 7 cm, (8) 

which is likely much larger than n c , and is thus 3> r co . A lower pulsar mass results in a higher 
ratio of r m /n c (Fig. 1). This means that the accretion with a rate of even Al ~ 5 x 10 15 g/s 
can not drive matter to reach the pulsar surface, but is in a propeller phase. That r m r\ c 
could be a direct cause for PSR B1259-63's being a radio pulsar. 

2 Even in this case of negligible momentum of accreting wind, a disk might form due to the interaction 
between the infalling medium and the rapidly spinning pulsar's magnetosphere, by which the rotational 
angular momentum of the pulsar is transferred to the accreted material. 

3 Note that Vfc(outside) ~ 98 km/s < u per ia < Vfc(inside) ~ 201 km/s in this case. 

4 This interaction transfers angular momentum from pulsar to accretion matter, and thus favors disk 
formation. 
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Fig. 1 The accretion rate derivative (solid lines, in unit of g • s -1 • P^b); is 
computed as a function of pulsar's mass in order for the observed braking index to 
be -36.7 by Ea. (|10|l . Each curve is for a different value of e, the meaning of which 
can be found in Eq.©. In the calculation, we choose an average accretion rate of 
M = (20 days/P or b)Mo, assuming the pulsar accretes for 20 days with a rate ~ Mq. 
The dotted line shows the magnetospheric radius r m in units of the light cylinder 
radius n c . It is shown that r m /r\ c ~ 5.5 for M p — M©, but increases as M p decreases. 



The magnetospheric interaction between a pulsar and its surrounding disk is still not well 
understood. Menou et al. (2001) discussed disk torque for the case of r m ~ r lc , and recognized 
that the torque would likely be largely reduced or entirely suppressed if r m ^> r\ c . Both the 
magnetodipole radiation and the propeller effect may result in the pulsar's spindown, 

n = -An 3 - eBM 3/7 n, (9) 

with 

A = 2/i 2 /(3c 3 /), B ~ 2.88 x 10 26 Mp/ 3 /Z, 

where I = 2M p R 2 /5 is the moment of inertia, and the propeller torque, Cl = — IMr^ljl of the 
second term in Eq.10, proposed by Menou et al. (1999) is applied, with a suppression factor 
introduced, e. The magnetodipole torque is generally smaller, but almost of the same order, 
as the propeller torque for e = 1. In other words, the magnetodipole torque is generally larger 
than the propeller torque for e < 0.1. The calculation of /x in Eq.Q is therefore consistent. One 
can also calculate the braking index from Eq.©, 

n = ™ = -^-[(3Ail 2 + eBM 3 / 7 )h + ^eBSlM-^M}. (10) 

Since 51 < 0, it is evident that the braking index n > for an accretion with constant rate, 
M = 0. However, the observed index is negative, n — —36.7 < 0, which means the effective 
accretion rate increases, M > 0. 

Since the accretion rate derivative, M, and the suppression factor, e, are two parameters 
hitherto unknown, let's find the relationship between them in order to explain the observed 
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braking index n — —36.7. The results, numerically calculated via Ea. (|10fl . are shown in Fig. 1. 
It is found that high e and/or low pulsar-mass (M p ) would result in a small M. The averaged 
accretion rate of M = (20 days/P orb )Mo ~8x 10 13 g/s. PSR B1259-63 could be a low-mass 
quark star (with M < O.1M ) if e < 10~ 4 and M < 10 12 g/s/P orb . 

The change of accretion rate on long timescales (<; P or b) was considered in the previous 
paragraph. It is natural to suggest that the accretion rate also increases before (but decays 
soon after) the periastron, e.g., exponentially, in a form of M ~ Moe ± */ T on a short timescale 
t < Porb- This short-term variation of M could significantly affect the timing behaviors close to 
the periastron, which could be the reason that a few timing points have systematic discrepancies 
at periastron passages. 

4 CONCLUSION AND DISCUSSION 

A few scenarios for understanding the x-jumps at the periastron are presented, including planet- 
like objects around the Be star, asymmetric mass-distribution by the Be star's oscillations, and 
a density-wave pattern in the Be star's disk excited by objects orbiting the Be star. The first 
one could be more likely than the others. We show that an accretion disk forms around the 
rapidly rotating, strongly magnetized pulsar because of the circularization radius being much 
larger than the radius of the pulsar's magnetospheric radius. The inner region of the disk is 
also beyond the light-cylinder, provided that the pulsar's accreted matter is replenished in the 
Keplerian disk of companion SS 2883. Propeller torque would therefore result in the pulsar's 
spindown, besides that due to magnetodipole radiation, though it may be suppressed by a factor 
e. The accretion rate should increase on a timescale of <; P or b in order to explain the negative 
braking index n = —36.7, but the reason for this increase is not certain (for example it may 
be due to a secular changing of a p via tidal effects) . Timing models with the inclusion of disk 
propeller effects would be necessary for fitting the timing data, and further observations would 
in turn check the model. 

The PSR J0045-7319 binary system (Kaspi et al. 1996) is similar to PSR B1259-63, in that 
tidal interactions may also be enhanced near the periastron in that system (Witte & Savonije 
1999). But why are their timing behaviors so different? This could be due to (1) the differences 
of pulsar masses (the mass of PSR J0045-7319 might be much lower than that of PSR B1259- 
63, M1259 <C M0045), (2) the differences of companion masses (Mf 259 > Mq 045 ), and (3) the 
different evolutionary stages of the companions (a younger Be star may oscillate with higher 
amplitude, or, planet-like objects could have been captured or evaporated by an old B star). 

A very negative braking index might be a feature of the disk with increasing mass. Long- 
term timing of AXP/SGRs and CCOs is proposed to see if there are propeller disks around 
them (the index n > for a decaying disk). In addition, this work has implications for studying 
long-term timing behavior of radio pulsars with propeller disks formed from either supernova 
fall-back onto young neutron stars or interstellar medium medium accretion onto old ones 
(probably near the death line in the P — P diagram) . 
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